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1 Introduction and main results 

We consider the spin 1/2 XXZ Heisenberg model on the dimensional lattice 
Z"^. For any finite volume A C Z*^, the Hamiltonian is given by 

Ha = - Yl ^^'(s^^s^^ + s^s^) + (1-1) 

\x-y\^l 

where A > 1 is the anisotropy. We refer to the next section for more precise 
definitions. By adding an appropriate boundary term one can insure that the 
ground states of this model describe an interface in the (1, 1, . . . , 1) direction 
between two domains with opposite magnetization. For a particular choice 
of boundary term, the model has exactly one ground state tpn for each fixed 
number of down spins, n. We call these the canonical ground states. In 
analogy with statistical mechanics of particle systems one can introduce the 
grand canonical ground states of the form 

n 

It turns out that these states are inhomogeneous product states 0]. In this 
paper, we consider a class of perturbations of these product states, of which 
we calculate the energy. By the variational principle this leads to bounds for 
the energy of the first excited state of the model. As the excitation spectrum 
above the interface states is gapless this bound should vanish as the 

volume tends to infinity. This is indeed the case (see (|1.2| )). 

The perturbations we consider are in correspondence with functions / : 
A ^ C. Furthermore, we consider functions which are slowly-varying in 
all directions perpendicular to (1, 1, . . . , 1) though they may have discrete 
jumps parallel to this direction. In other words || V/ • w||oo -C ||/||oo for all 
t>±(l,l,...,l). We consider general perturbations of this type and conclude 
that the optimal perturbations, in the sense of minimizing energy, are local- 
ized near the interface. With this restriction, the Hamiltonian, projected to 
and restricted to the appropriate subspace, is just the Laplacian 

This result may be compared to the recent bound of 0. The main dif- 
ference is that there we considered a canonical ensemble, for which there 
were a fixed number of down-spins (hence a fixed number of up-spins). We 
developed a version of equivalence of ensembles whereby we estimated the 
canonical expectation of a gauge invariant observable by a grand canonical 



3 



expectation, provided that the interfaces of the canonical and grand canonical 
states occupied the same position. 

In the present paper, we begin with the grand canonical ensemble, so that 
we make no reference to equivalence of ensembles. Specifically, we consider 
a cylindrical region of total height L + 1 and whose cross-section is a region 
Qji with linear size R. Then a class of excitations is parametrized by smooth 
functions $ on a fixed domain Q = R~^Qii. 

Main Result: Excitations on A have a normalized energy 

{¥\H\^f) 1 w^ni^n) , 



lL2(n) 



where 

^L/2-l 



^(A,/i) 



^«=-L/2 sech(a[/ — /i]) sech(Q;[/ + 1 — /i]) 
Yl^=-L/2 sech(a[/ — /i]) sech(a[/ — /x]) 



-L/2 

Here, n is a real parameter of the grand canonical ground state describing 
the location of the interface between the regions of homogeneous up and down 
spins. As fi ^ —oo, the ground state has all spins up, and for fi oo, all 
spins are down. For all /i G M, and sufficiently large L, g satisfies the bounds 

1 



— <A-v/A^<(7(A,^)<l 



Remark: The normalized energy of ( |1.2|) is exactly the same as that 
for the Laplacian. Equating the first variation to zero, we see that the local 
extrema of the normalized energy are precisely the solutions of = — A$ 
(here is the Laplacian), and A = || V$||^2(Q)/||$||^2(f^)- The space of 
excitations we consider does not form an invariant subspace of H, so that 
the eigenvectors of the Laplacian are not truly eigenvectors of H. But, using 
the variational inequahty, we see that the spectral gap of H is bounded thus: 

7i<^-9(A.,)(l + 0(i,)), 



where Ai is the first positive eigenvalue of — with Dirichlet boundary 
conditions on the domain Q. 
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2 The Spin-| Heisenberg XXZ Ferromagnet 

A quantum spin model, such as the Heisenberg XXZ ferromagnet, is defined 
in terms of a family of local Hamiltonians H\, acting as self-adjoint linear 
operators on a Hilbert space TYa. This family is parametrized by finite subsets 
A c Z'^. 

We choose A to be "cylindrical" in the following sense: Let {e^}^^]^ be 

the set of coordinate unit vectors and define the vector = Yl'^j=i ~ 
which is the axial direction for the cylinder. Define the func- 
tional l{x) = X' ■ = Yl'^=i^''^ where x = Yl'j=i^''^j- Observe that the 
kernel of I in is a (rf — 1) -dimensional sublattice perpendicular to the axial 
direction. Take for the base of A a finite subset of this {d — l)-dimensional 
sublattice, and call it F. A discrete approximation to the line of all scalar 
mutliples of is the one- dimensional stick S. S is a bi-infinite sequence of 
points {xn}'^=-oo such that xq = and all other points specified by 

the relation x^ - x^-i = e^modd- So 

E = {. . . ,-{ed + Cd-i^ \-ei + ed),-{ed + ed-i^ h ei), . . . , -e^, 

0, ei, {ei + 62), . . . ,{ei + e2-\ h e^), (ei 62 H \-ed + ei),.. .}. 

A finite stick of length L + 1, where L is even, is T,l — {x E T, : —L/2 < 
l{x) < L/2}. Now define A to be the translates of F along E^, i.e. 

A = F + Ei = {a;-H2/:a;eF,|/eEi}. (2.1) 

Let us now define nearest neighbors to be points x,y G Z,'^ such that 
|/(a;) — l{y)\ = 1 and \\x — y\\ii = 1. Also, we define oriented bonds between 
nearest neighbors as ordered pairs {x,y) satisfying l{y) = l{x) + 1 and \\x — 
y\\ii — 1. Hence {{x,x + ej)}j^i is the set of all oriented bonds with lower 
point X. The collection of all oriented bonds with both points in A, will be 
called B{A). 

The local Hilbert spaces are TYa = (C^)®'"^'. Each copy of comes with 
an ordered basis (|t) , and a spin-i representation of SU{2) defined by 
the Paufi matrices: 

~VV2 yl' "VV2 yl' ^ \ -1/2 J ■ 

(2.2) 
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(So, for example, S^^^ |t) = | |t) and 5*^^^ |J.) = — | ||).) We consider a family 
of Hamiltonians parametrized by a real number A > 1. In order to define the 
total Hamiltonian, we first define local Hamiltonians h^y for each oriented 
bond (x, y): 

h.,y = -A-^(^«4^) + 4^)) - ^i^)^^) + \ + ^A(A)(4^) - ^(^)), (2.3) 
where ^(A) = |a/1 — 1/A^. The total Hamiltonian is 

Ha= ^Ir (2.4) 

(a;,y)GB(A) 

A parametrizes "anisotropic coupling". The case A = 1 is the isotropic 
model, also known as the Heisenberg XXX ferromagnet, which exhibits SU (2) 
symmetry (because H\ commutes with 5"^, 5"^ and S^). 

We find it convenient to introduce a positive constant a, which solves A = 
cosh(a). We note that the nearest neighbor interaction h^y is an orthogonal 
projection 

hxy = \^xy) {ixy \ ® ]lA\(x,y), (2.5) 

where 

. _ e-"/^|iT)-e"/^|n) 

^xy - r- ■ K^-'O) 

y 2 cosh(aj 

This also shows that each hn^y is a nonnegative self-adjoint operator, hence 
i^A is, as well. To simplify the notation we will often drop the subscript A 
when the volume is obvious from the context. 



3 Ground States and a Perturbation 

The ground states of the XXZ ferromagnet can be calculated exactly [Q. We 
will choose a particular ground state and construct an orthogonal subspace 
(but not the entire orthogonal complement) which is parametrized by H^- 
functions on a compact domain f2o C M*^^^. The inner product becomes 
approximately the inner-product and the orthogonal projection of the 
Hamiltonian is approximately the Laplacian. 

The lowest eigenvalue for if, which is zero, has a (| A| + l)-fold degeneracy 
in the eigenspace. This space of ground states is spanned by the simple tensor 
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ground states, which we will call grand canonical states. Specifically, let z 
be any complex number, and /i = Re(2). Define the vector 

Vx\z) = ^ , (3.1) 

a/2 cosh(a[Z^ - jj]) 

for each site x G A. We define the product of these vectors 



Mz)=Q^Vx{z), (3.2) 

and we may quickly establish that it is a ground state. Indeed, the oriented 
bonds are defined between points x and y with l{y) = l{x) + 1, from which 
we see 

(Ti Kiz) Vyiz)) = e"(iT \vx{z) ® Vy{z)). (3.3) 

This implies Vx{z) ®Vy{z) is orthogonal to ^xy, for each {x^y) e -B(A), which 
proves that 'ipQ^z) is a ground state. As we have said, the states iIjo{z) span the 
entire ground state space, as z ranges over all the complex numbers ||^. (More 
than this can be said. The simple tensor ground states are parametrized by 
elements of CP^, so that the submanifold of all such states in Ti is topolog- 
ically a sphere. But to obtain the north and south poles of the sphere, it is 
necessary to take the limits z — > cxd and z — *• —oo.) 

Let us now fix z, and for simplicity we will just write ipo and Vx without 
explicit reference to z. For each site x we define a vector orthogonal to Vx-, 

Wx = ^ • (3.4) 

^J2 cosh(a[/^ - /i]) 

We will make use of Wx to define an orthonormal system of states 



V=Wx® Vy, (3.5) 

y^A\x 

where x ranges over A. Each of these states is also orthogonal to ipQ, let 
us call their span V . An arbitrary state in V is characterized by a function 
/ : A ^ C. Exphcitly, ^Ij^ = T.x&Kf{^W- It is then clear that (^/'^|^/'^) = 
T.x<^Kfi.^)9{x). 

Our interest is the case that A y Z"^, i.e. the thermodynamic limit. 
In terms of Vx and Wx-, we see that the local interaction h^y describes a 
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nearest-neighbor interaction. It may be interpreted as a bilinear form, which 
is a first order finite-difference operator in each variable. To be clear, a 
straightforward calculation gives 

1 

\i' {hxylip^) — - sech(Q;) sech(ci;[Z3. — //]) sech(Q;[Zy — //]) 

x(cosh(a[/j^ - fi])f{y) - cosh(a[/^ - 

x(cosh(a[/j^ - l4)g{y) - cosh(Q;[/3; - i^])g{x)). (3.6) 

Recall that /i = Re{z)) and the energy is 

L/2-1 d 

{^pf\H \r) = Yl E E (^^1 1^^) , (3.7) 

l=-L/2 xeVi 3=1 

where F; refers to the set of points x e A with l[x) = I. In the thermodynamic 
hmit, we may scale the plane — {v & : v-e^ — so that H becomes, to 
first order, a differential operator with respect to each direction of the plane. 
However, the inhomogeneity in the e^, direction admits no such scaling for 
that coordinate, so that H is genuinely a finite-difference operator even in 
the thermodynamic limit. 

This intuitive description of the last paragraph is made precise, now. 
Let n be a bounded, open subset of with a boundary. Let be 
the dilation R ■ Vt = {Rx : x E Q}, and let F = fl Z*^ be the discrete 
approximation to Qr. As before, F is the base of A. Now we choose a smooth, 
complex- valued function $ on Q, and extend it to the infinite cylinder fix Re* 
so that V$ • e* = 0. (In other words, $ is constant along the direction 
e^,.) Let (f){x) — ^{x/R), which is defined on Qji x Re* with the property 
that V0 • e* = 0. Finally, let f{x) = F{l,j.)(j){x)^ where F is a sequence 
F{—L/2), . . . , F{L/2). Note that / is not the most general form possible for 
a function on A, most notably because it is the product of functions which 
vary on perpendicular subspaces. However, the span of such functions does 
correspond to all of V for a fixed value of L and R. 

Next we consider the norm and energy for such a state. We will introduce 
estimates for these quantities, but we will postpone the actual error terms 
until the next section. First we replace the sum over F with the integral over 
fi, and thus obtain an expression for the norm: 

L/2 

ii^'i^') = E E 1/(^)1' 

l=-L/2 xeVi 
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LI1 



- iri Ej^(or-;;^/j*(-)r^- (3-8) 

To obtain an approximation for H j'?/'-^), we decompose a step of / along 
a coordinate direction into a step parallel to e^, and a step perpendicular to 

/(x + ej) = + + e,) 

^ + l)0(x) + + l)V0(x) ■ e,. 

Then using the fact that 

d 

^ V0(x) ■ Cj = V0(x) ■ = 0, 
i=i 

and referring to ( |3.6| ) and (|3.7|) we have the apparently cumbersome expres- 
sion 

I 1^ / 2cosh(a) m{Q) ^ ' ^ ^' 

L/2 

X |^sech(a[/ — /i]) sech(a[/ + 1 — /i]) 

Z=-L/2 

I cosh(a[/ + 1 - + 1) - cosh(a[/ - fx])F{l)\^ 



2it^cosh(a) m[il) Jq 



L/2 

X ^ sech(a[/ - /i]) cosh(a[/ + 1 - + 1) 

Z=-L/2 



We notice that the first summand is order 1, while the second summand 
is order 1/B?. We wish to minimize the energy in the limit — > oo, so it 
seems sensible to eliminate the order 1 summand. This is accomplished by 
letting F{1) = ^ sech(a[/ — /i]), or any constant multiple thereof. One point 
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of interest is that the perturbation takes place primarily in a neighborhood 
of the interface. The expression for the energy is 

L/2-1 

X sech(a[/ — yu]) sech(a[/ + 1 — yu]). (3.9) 

l=-L/2 

Similarly, (|3.8| ) may be rewritten as 

in 1 r 

4 m(fi)7^ ^^^^^ 

Taking the ratio, we arrive at a normalized energy 
{ipf\H\ij^) sech(a) l|V$|| 



|2 



Yld=-Li2 sech(a[/ — /x]) sech(a[Z + 1 — /x]) 

X ■ 

X]/=-L/2 sech(a[/ — /i]) sech(a[/ — yu]) 



(3.11) 



Let P be the orthogonal projection to the subspace of perturbations con- 
sidered so far, i.e. the span of , where f{x) = |sech(a[/a; — /i])0(x). Then 
the projection of H to this subspace is PHP. We have determined that 
PHPipf = ip^ where g has in place of $ 

sech(a) sech(a[/ - /i]) sech(a[/ + 1 - /i]) 

2-R X]/=-L/2 sech(Q;[/ — yu]) sech(a[/ — yu]) 

(We write for the Laplacian. The symbol A is reserved for the anisotropy.) 
We should note that it really is necessary to consider PHP instead of H. 
The reason for this is that 

^ —2 cosh(a[/^ — fj])wx <^ Vy + 2 cosh(a[/j^ — fi])v,j. ®Wy + 2 sinh(Q;)w^ ® Wy 
^ a/2 cosh(a;[/2,. — yu]) ■ 2 cosh(a[/y — yu]) ■ 2 cosh(a) 

(3.13) 

which means that H does not preserve the total number of f^'s or Wx^s. Thus 
the perturbations we have considered (those with a single Wx) do not form 
an invariant subspace of H. 
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4 Error Terms 

We now come to the task of tying-up some loose ends, in order that non- 
rigorous approximations can be replaced by rigorous bounds. We start with 
a simple lemma. 

Lemma 4.1 Let F be a finite subset of a lattice L. Let Q be the Voronoi 
domain of V with respect to L, and let be the Voronoi domain for the 
single site e L. Then, for a smooth function : Q — > C, 

[ 4>{x)dx <\\d^4>\\op,oo-—^ I ^o?x, (4.14) 

where is the second-derivative matrix and 

||a^0||op,oo = sup sup (4.15) 



Note that the second moment m(i7o) ""^ / l^cprfa; is bounded by the radius of 
the Voronoi domain, which is in turn bounded by the distance of nearest 
neighbors of L. 

Proof: For the Voronoi domain Qq of 0, we observe that 
1 r , . , 1 



[ 0(x)dx-0(O) = -^TTT [ [Hx)-m]dx 



1 r 

V0(te) ■ xdtdx 



0) jQo Jo 
1 r rl rt 



m(f2o 



X ■ V (j){sx)x ds dt dx 



0) Jno Jo Jo 

+— ^ / ^(p{0)-xdx 
m{no) Jno 



— T— ^ / / (1 — s)x ■ d'^(j)(sx)x ds dx 
m{^o) Jno Jo 

+V0(O)-— ^ [ xdx. 
m{flo) Jn^ 



But the centroid of is 0. Thus 



1 



^(^o) Jno ^(^o) Jno 2 



1 f \x 



|2 



(j){x)dx-(l){0) <—— h^dx X \\d'(l)\\op,oo. (4.16) 
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The lemma follows by decomposing Q into the |r| affine copies of Qq, one for 
each site, and adding the inequalities obtained from ([4.16|) . g 

Using the result of this lemma, we make rigorous the approximation of 



dp). Thus 



where 



hi < ^1191* 



I op,oo- 



(4.18) 



(We have used the fact that the distance between nearest-neighbors for F is 
■\/2.) In order to fix the approximation of p.9| ), we begin with the elementary 
bound |0(x + Cj) — 0(x) — V0(x) • ej\ < \\\d'^(f>\\op,oo and its natural successor 



I J2 I'Pi^ + ^j) - - l|V0(x)||2| < f/(||V0||oo + ^\\d'nop,oo) \\d\ 

Using this estimate, as well as the lemma, we may replace ( |3.9|) with 

11 L/2 



^11 op,oo- 

(4.19) 



SR^cosHa) ^ sech(a[Z-/.])sec(a[/ + l-/.]) 

^ ^ l=-L/2 



where 



and 



e.|<|(l|V*IU + j^l|#t 



f3l<;jpll#|Vt| 



I op,oo 



I op, oo ■ 



(4.20) 

(4.21) 
(4.22) 
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